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We give a comparison of the spectrum of Yang-Mills theory in D = 3 + 1, recently derived with 
a strong coupling expansion, with lattice data. We verify excellent agreement also for 2 ++ glueball. 
A deep analogy with the D = 2 + 1 case is obtained and a full quantum theory of this approach is 
also given. 
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Yang-Mills quantum theory in the low energy limit represents one of the most difficult physical problems to face 
with. Difficulties in managing this theory have largely motivated the use of large computers to get the spectrum. 
Recently, after a relevant theoretical work mostly due to Nair and Karabali [l], 0, H, 0, [B| , it has been possible to 
obtain the spectrum of this theory in 2+1 dimensions [f|. The agreement is quite satisfactory and improvements on 
this line are possible as further clarifications on the lattice data are put forward. This approach seems promising also 
to manage the 3+1 case p}- 

Recently, we proposed an approach to treat the 3+1 case by a stron g co upling expansion, then we showed that this 
expansion is nothing else than a gradient expansion 0, 0, [l(J, 11 , 12l Il3l . This approach gives the spectrum of the 
++ glucballs in good agreement with recent lattice computations |l4 , Il5j | . Furthermore we were able to give a very 
nice agreement between our propagator [l(| and lattice data [l6| when the gluon mass is taken to be 389 MeV that 
is the only parameter to fit. 

In this paper we have two aims. Firstly a full quantum theory is give n using functional integration. Secondly, we 
verify that the spectrum of the glueballs obtained on the lattice [Til 15 1 is properly reproduced. As a by-product we 
obtain a deep similarity with the aforementioned D=2+l case. 

In order to make this paper self-contained we give here a derivation of the propagator. The starting point is a A0 4 
theory with the generating functional 

Zb1 = /[#]e{'-/ rfD 4i^^) 2 -^ 4 ^]} e {- i/dD 4* W)2+ ^ 2 ]} (1) 

where we have properly separated the spatial part from the time derivatives of the field. A gradient expansion means 
that we make an expansion as 



Z[j) = cxp 



being 



z [j] = f [d0] e {^ d£ Ml(^) 2 -^ 4 +^]}. 



Zo\j] 



(2) 
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The final point is to write this latter integral in a Gaussian form. In order to realize this point we use the fact that a 
small time expansion exists for the solution of the equation 



J 



(4) 



after the change to dimensionless variables is made through /ij as i -» (J-o x , 2 ~~ * Mo D( f >2 an d introducing the 
4_jj . The small time expansion permits us to write the solution of this equation by a Green 

Mo 



coupling constant g 
function given by 



G + gG 3 = S(t) 



(5) 



as 



17, 1 



dt'G{t-t')j(t') + a / dt'G(t-t')(t-t')j(t') + b / dt'G{t-t'){t-t') 2 j{t') 
Jo Jo 



(6) 
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being a, b, . . . coefficients to be computed and dependent on the initial value j(0) and its derivatives. We recognize 
at the leading order the familiar result of Green functions for linear systems and this implies that, as a leading 
approximation, one can write a Gaussian expression for Z [j] 0. In order to write this expression, let us write down 
the Green function for our theory at the leading order. One has 



G{t) = 9(t) 



sn 



and its time reversed G(—t). Here one has sn a Jacobi elliptical function and 0(t) the Heaviside function. So, one can 
introduce the Feynman propagator A(x, t) = 5 D ~ 1 (x)[G(t) + G{— t)] and write down the Gaussian expression Q 



z o\j] = exp 



d xid x 2 j(xi)A(x 1 - x 2 )j(x 2 ) 



(8) 



and then the theory can be immediately integrated giving the two-point function (undoing dimensionless changes) [l 

B n 



A(p) 



E 

71 = 
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mi - Mo + «e 



being 



B n = (2n + 

and the mass spectrum of the theory given by 

m„ = (2n 



7T 2 (-l) n + 1 e -("+3) 7 ' 
] K 2 {i) l + e -(2n+l)7r 



1) 



A 



2K(i) \ 2f 4- D 

proper to a harmonic oscillator. In this case, the mass gap in the limit A 
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(9) 



(10) 



(11) 



(12) 



corresponding to the choice n = 0, that is produced by the self- interaction of the scalar field. In order to get eq.® 
we have used the known relation for the Jacobi elliptical function (lil ] 



sn(u, i) 



(— i)™ e -(™+3) 7r 



2tt 



K(i) ^ l + e -(2n+l) 7 r 
v ; n— 



sm 



(2n + l)- 



(13) 



do 



being K{i) = f » -~ 

u Y l+sin^ y 

The interesting case, D = 4, gives 



2if(i) 

1.3111028777 a constant. So, the spectrum is defined through a renormalized mass 8s 



2K(i) 



Mo- 



(14) 



At this order the theory appears finite but it is possible that higher order corrections to the Gaussian approximation 
could introduce the need for renormalization in this limit. Presently, we content ourselves with these results and we 
show how they apply to a Yang-Mills theory in 3+1 dimensions. 

A spectrum for the Yang-Mills theory can be obtained if we are able to map its quantum theory on that of the 
scalar field discussed above. This in turn means that we have to prove that a gradient expansion, identical to the 
one of the scalar field, does hold also for a Yang-Mills theory. For our aims, we turn to the classical theory. If we are 
able to prove the existence of the correspondence of the two theories in this case, the existence of the corresponding 
gradient expansion for the Yang-Mills theory is obtained. This we proved in Indeed, the classical scalar theory 
has the following Hamilton equations working now in D = 3 + 1 



d t <t> 
d t ir 



(15) 



V 2 0-Mo^-^ 3 
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We do the following changes and expansions 

= VAi (16) 



r 



tt = VA f TrW + W 1} + tW (2) 
\ A A^ 

4> = (o) + V 1} + i^ (2) + --- 

A A 

and then we get the following non-trivial equations for our expansions 



d T ^ = 7T<°> (17) 

d T <f>& = ^ 

9 r7r (°) = -0(°) 3 



that we recognize as a gradient expansion as should be. These equations can be solved and give a meaningful gradient 
expansion for a classical A</> 4 theory. We can proceed in a similar way for a SU(N) Yang-Mills field. The Hamilton 
equations are, in the gauge Aq — that needs to be fixed to carry on our computation, 2(J 21 1 



dtA a k = F§ k (18) 
%F« k = diF^+gf^A^ 



being g the coupling constant, f abc the structure constants of the gauge group, Ffa = diA k — dkA" + gj abc A\A c k and 
the constraint e^i 7 ^ + gf abc A k FQ k = does hold. So, let us introduce the following equations, as done for the scalar 
field, 



t = gt (19) 

ia(2 
Ok 



^ofe - 9*ok + t ok + ~*i 
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p a _ p a(0) 1 p o(l) , 1 w a(2) 
*lk ~ b lk + ~*lk + ~if lk 

At = Af ] + - g Af ) + ±Af } - 



Then, one has the perturbation equations 



d T Af> = F T } (20) 

1 a (!) _ Tpa(l) 



d T A^> = F { 



OA- 



g T P*(°) = jabc jcde ^d(0) ^(0) 

d F a< -^ = f abc f cde v4 b(1) Af^ ' A e ^ 0) + f abc f cde A b( - ^ ' A d ^ A e{ -^ + f abc f cde A h ^ A d{ - 0) A e ^ 



-f abc di (, 



^(0) A c(0)\ + fabcA b(0) /Q^ciO) _ Q^ciO) 



and we can recognize at the leadin g or der the homogeneous Yang-Mills equations. These equations display a rich 
dynamics as e.g. Hamiltonian chaos 2^, 23, 24 1. But chaotic solutions do not give a meaningful quantum field theory. 
But we easily realize that taking all the components of the Yang-Mills field as equal [2l| we can make a complete 
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correspondence between the gradient expansion of the classical scalar field (JT7J) putting A = Ng 2 , the so called t'Hooft 
coupling. An arbitrary mass is needed in this case because no mass scale exists in the Yang-Mills theory. This is 
an integration constant in this theory and is generally obtained, also on the lattice by fixing the string strength, by 
comparison with experimental data. In the following we will fix this constant to 1J, [l5( 



Ng 2 



A = 440 McV. 



(21) 



This choice, as we will se below, permits us to reproduce exactly the lattice spectrum. So, as this series exists for the 
scalar field it does exist for the Yang-Mills field [9(] . We used the fact the g in D = 3 + 1 is dimensionless as happens 
to A for the scalar field and we have used it as a useful placeholder for our expansion. But the same result could be 
obtained by removing it from the theory using the duality principle in perturbation theory [25| . 

Now, in order to complete our proof, we consider this procedure starting from the functional integral of the Yang- 
Mills theory. One has, choosing now the gauge d^A a ^ = 0, [26j 



being 



and 



Zym[j] = / [dA][dc][dc]e 



S = / d 4 



iSa+iSiNi 



-d li A®d fJ 'A au + d^c a d^c a 



S 



INT 



= d 



gf abc d ll AZA btt A a ' + a —f hc f ars "A\A%A rii A s 



gf abc d^A hil c c + j atl A a u 



(22) 



(23) 



(24) 



We can assume for A't that one of the spatial components is zero and the others are all equal in agreement with the 
classical case. This is not a truncated Yang-Mills theory but a proof of the existence of a quantum solution completely 
mapping a Xcj) 4 theory that holds in the infrared. For consistency reasons we do the same for the ghost field that in 
this way happens to decouple from the Yang-Mills field and obtaining in this way a propagator going like 1/p 2 that 
is diverging in agreement with lattice computations [3]. So, finally one has 



S = -(A 2 -1) J 



d A x 



2 M 



dede 



and 



Sint = —(A 2 - 1) 



d 4 x 



Ng 2 



A + j A 



(25) 



(26) 



where we note the factor A 2 — 1 needed to preserve the number of degrees of freedom. Then, a gradient expansion 
holds also for the quantum case and the given theory exists being the same as for a X<f> quantum field theory in 
the infrared. Having a zero mass theory, scale invariance is retained as for the Yang-Mills theory. Then, as for the 
classical case, we need an arbitrary mass scale to normalize the theory. In this case, as noted above, we call this mass 
A. 

The above correspondence in a gradient expansion between a scalar field theory in D = 3 + 1 and a Yang-Mills 
theory gives us all the information we need about the spectrum and the propagator of the latter. Indeed we have 



A 



ym(p) = Y] — 



B n 



n=0 J 



M(n) 



being 

B n = (2n 

and the mass spectrum of the theory given by 



TT 2 (-l)"+ 1 e -(«+3) 7r 
] K 2 (i) 1 +e -(2n+l)7r ' 



M(n) 



(2n + l) 



2K(i)' 



(27) 



(28) 



(29) 
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Excitation 


Lattice 


Theoretical 


Error 


a 




1.198140235 




0++ 


3.55(7) 


3.594420705 


1% 


++ * 


5.69(10) 


5.990701175 


5% 


0++" 




8.386981645 





TABLE I: Comparison for the ++ glueball 



Excitation 


Lattice 


Theoretical 


Error 


* 

a 




2.396280470 




2++ 


4.78(9) 


4.792560940 


0.2% 


2++* 




7.188841410 





TABLE II: Comparison for the 2++ glueball 



The mass gap is 



M(0) - M®*' m 



to be identified with the gluon mass dynamically acquired by the strong interaction due to the self-interaction terms 
of the Yang-Mills field. We note that, differently from the case D = 2 + 1, the string strength is an arbitrary constant 
that should be derived by experiment or a higher order theory. Indeed, we are not able to apply this correspondence 
to a lower dimensional Yang-Mills theory, this is indeed a lucky case, and then the approach given in [J, [3, 0, 0, HI, 0] 
is welcome as could be its application to a D = 3 + 1 case to be compared with our results. 

We now compare our results with lattice data 14, TH(3- For the propagator we just note that the agreement is 



excellent fixing the gluon mass to 389 MeV [lfjj. We do not pursuit this matter further here waiting for other data 
from lattice computations. We note that we have obtained a result that does not depend on N, the group order. This 



is aligned with lattice results as given in 14j,[l5|. So, for our convenience we limit the comparison to the continuum 
for the SU(3) group. Lattice data in L4J, ll5| are given only for ++ , ++ * and 2 ++ and we assume these results are the 
ones to be more confident. Further results for the spectrum can be read in [27| and (28|, the latter being an unofficial 
extension of the results given in 14, 15| . We do not aim to go that far but just to have a proper understanding of 



how to get, given the propagator, higher order excitations. Initially we just notice that M{n) is indeed the spectrum 
of ++ . The reason for this relies on the fact that these are poles of the propagator and represent physical states of 
the theory. So, we have to interpret them in this way and indeed we obtain agreement with lattice computations. In 
Tab. |T]we can see an exceedingly good agreement between lattice and theoretical data. The interesting point is about 
the excitation called a in our table. We called this particle in this way as M(0) = 527 MeV that is about the mass 
of the a resonance in (29l | . This is another theoretical result of our approach and currently it is not expected to be 
seen in lattice computations. We expect to see its excited states in higher excited states of the spectrum. 

Our aim and the main motivation of this paper is to see what happens when higher excited states than ++ arc 
considered. This analysis can be carried out as already done in D=2+l [||. We consider a time independent correlator 
given from the propagator (|27p by taking pq = and Fourier transforming in space obtaining 

+ oo . 

that has an identical form to the asymptotic correlator obtained in D=2+l for ++ states as should be[(|. Taking 
the square of this expression we obtain for the 2 ++ spectrum 

M ^ m) = MM±MM = (2n + 2m + 2) * (32) 

V " V a 2K(i) 

that gives Tab [LT1 through M2(0, 1) = M%(1, 0). We observe a degeneracy, the same of the two-dimensional harmonic 
oscillator. The agreement is again astonishingly good. 

A few words should be spent about the strong similarities between the D=2+l and D=3+l cases. In both cases we 
observe the same asymptotic form of the propagator and the same harmonic oscillator spectrum that gives exceedingly 
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good agreement with lattice data. As a final note we point out that some analysis on the lattice support the view 
that the spectrum of the theory in D = 2 + 1 is indeed that of a harmonic oscillator 3(J 31 1 . 



We have given a full account of a strong coupling approach that gives the spectrum of a Yang-Mills theory in D=3+l. 
The agreement with lattice data are exceedingly good and the analogy with the Yang-Mills theory in D=2+l is really 
striking putting our results on a sound ground. Finally we have also given a fully quantum field formulation of 
Yang-Mills theory in our approach. 
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